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ABSTRACT 

We have used the Very Large Array (VLA) to search for radio emission from 
the globular cluster Gl (Mayall-II) in M31. Gl has been reported by Geb- 
hardt et al. to contain an intermediate-mass black hole (IMBH) with a mass of 
~ 2 X 10^ Mq. Radio emission was detected within an arcsecond of the cluster 
center with an 8.4 GHz power of 2 x 10^^ W Hz~^. The radio/X-ray ratio of Gl 
is a few hundred times higher than that expected for a high-mass X-ray binary in 
the cluster center, but is consistent with the expected value for accretion onto an 
IMBH with the reported mass. A pulsar wind nebula is also a possible candidate 
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for the radio and X-ray emission from Gl; future high-sensitivity VLBI observa- 
tions might distinguish between this possibihty and an IMBH. If the radio source 
is an IMBH, and similar accretion and outflow processes occur for hypothesized 
~ 1000 Mq black holes in Milky Way globular clusters, they are within reach of 
the current VLA and should be detectable easily by the Expanded VLA when it 
comes on line in 2010. 

Subject headings: accretion, accretion disks — galaxies: individual (M31) — 
globular clusters: individual (Gl) — radio continuum: galaxies 



1. Introduction 



Most galaxies with massive spheroidal components appear to harbor central black holes 
(BHs), with masses ranging from a few xlO^ Mq to over 10^ Mq. These BH masses are 
well correlated with both the luminosity and th e velocity dispersion of the galaxy spheroid 
( iGebhardt et al.ll2000l : iFerrarese fc Merrittll2000l ). implying that the formation of the central 
BHs is connected intimately with the development of the galaxy bulges. However, a bulge 
is not necessarily a prerequisite for a massive BH. On the one hand, neither the late-type 
spiral galaxy M33 (Gebhardt et al. 2001) nor the dwarf spheroidal galaxy NGC 205 (Valluri 
et al. 2005) shows dynamical evidence for a massive BH. On the other hand, both the Sd 
spiral galaxy NGC 4395 and the dwarf spheroidal galaxy POX 52 contain active BHs with 
masses ~ 10^ Mq, although neither contains a classic bulge (Fillipenko & Ho 2003; Barth 
et al. 2004). Greene & Ho (2004, 2007) find that optically active intermediate- mass BHs 
(IMBHs, Mbh < 10^ Mq), while rare, do exist in dwarf galaxies, but optical searches are 
heavily biased toward sources accreting at high Eddington rates. Alternate search techniques 
are needed to probe the full demographics of IMBHs. 

Direct dynamical detection of IMBHs is currently impossible outside of the Local Group. 
However, Gebhardt, Rich, & Ho (2002, 2005) have found dynamical evidence for an excess 
dark mass of (1.8±0.5) x 10^ Mq a t the center of the globul ar cluster Gl in M31; the evidence 
for this IMBH was questioned by iBaumgardt et al.l (120031 ) . but supported by the improved 
data and analysis of Gebhardt et al. (2005). The physical nature of the central dark object 
is difficult to prove: it could be either an IMBH or a cluster of stellar remnants. While 
the putative presence of a BH in a globular cluster center may appear unrelated to galaxy 
bulges, the properties of Gl, including its large mass, high degree of rotational support, and 
multi-a ged stellar populatio ns, all suggest that Gl is actually the nucleus of a stripped dwarf 
galaxy (IMeylan et al.ll200ll ). Most intriguingly, the inferred BH mass for Gl is about 0.1% 
of the total mass, consistent with the relation seen for higher mass BHs, and consistent with 
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predictio ns based on mergers of BHs (IMiller &: HamiltonI l2002l ) or stellar mergers in dense 
clusters JPorteries Zwart fc McMillan|[2002 ). 



Recently, iPooley fc RappaportI (120061 ) reported an X-ray detection of Gl, with a 0.2- 
10 keV luminosity of Lx ~ 2 x 10^^ ergs s~^. Although this may represent accretion onto 
a central BH, it is within the range expected for either accretion onto an IMBH or for a 
ma ssive X-ray b inary. Unfortunately, the most accurate X-ray position determined recently 
by iKongI (120071 ) does not have sufficient accuracy to determine whether the X-ray source 
is located within the central core of Gl, which would help distinguish between these two 
possibihties. The radio/X-ray ratio of Gl prov i des a n additional test of the nature of the 
Gl X-ray source. As pointed out by iMaccarond (120041 ) and Maccarone, Fender, & Tzioumis 
(2005), deep radio searches may be a very effective way to detect IMBHs in globular clusters 
and related objects, since, for a given X-ray luminosity, stellar mass BHs produce far less 
radio luminosity than supermassive BHs. The relation between BH mass, and X-ray and 
radio luminosity empirically appears to follow a "fundamental plane," in which the ratio of 
radio to X-ray luminosity increases as the ~ 0.8 power of the BH mass (Merloni, Heinz, & di 
Matteo 2003; Falcke, Kording, & Markoff 2004). For an IMBH mass of 1.8 x 10^ M© in Gl, 
one thus would expect a radio/X-ray ratio about 400 times higher than for a 10 Mq stellar 
BH. In this paper, we report a deep Very Large Array (VLA) integration on Gl and a radio 
detection that apparently confirms the presence of an IMBH whose mass is consistent with 
that found by Gebhardt et al. (2002, 2005). 



2. Observations and Imaging 

We obtained a 20-hr observation of Gl using the VLA in its C configuration (maximum 
baseline length of 3.5 km) at 8.46 GHz. The observation was split into two 10-hr sessions, one 
each on 2006 November 24/25 and 2006 November 25/26. Each day's observation consisted 
of repeated cycles of 1.4 minutes observation on the local phase calibrator J0038+4137 and 
6 minutes observation on the target source Gl. In addition, each day contained two short 
o bservations of 3C48 (J0137-I-3309) that were used to calibrate the fiux density scale to that 



of iBaars et al.l (119771 ). Thus, the total integration time on Gl was 14.1 hr. We also obtained 
a total of 9.5 hr of observing in C configuration at 4.86 GHz on 2007 January 13/14 and 
2007 January 14/15, using a similar observing strategy, and achieving a total of 7.3 hr of 
integration on source. 

All data c alibration was carried out in NRAO's Astronomical Image Processing System 



(lGreisenll2003l ). Absolute antenna gains were determined by the 3C 48 observations, then 



transferred to J0038+4137, which was found to have respective fiux densities of 0.52 mJy and 
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0.53 mJy at 8.4 and 4.9 GHz. In turn, J0038+4137 was used to calibrate the interferometer 
amplitudes and phases for the target source, Gl. Erroneous data were flagged by using 
consistency of the gain solutions as a guide and by discarding outlying amplitude points. 

The VLA presently is being replaced gradually by the Expanded VLA (EVLA), which 
includes complete replacement of virtually all the electronic systems on the telescopes. Since 
antennas are refurbished one at a time, the VLA at the time of our observations consisted 
of 18-20 "old" VLA antennas and 6 "new" (actually, refurbished) EVLA antennas, having 
completely different electronics systems. Although all antennas were cross-correlated for our 
observations, we found subtle errors in some of the EVLA data. Thus, to be conservative, 
we discarded the data from all EVLA antennas except for 3 antennas that were confirmed 
to work very well on 2006 November 24/25. 

The radio data were Fourier transformed and total-intensity images were produced in 
each band, covering areas of 17'xl7' at each frequency^] These images were CLEANed in 
order to produce the final images. At 8.4 GHz, the rms noise was 6.2 /iJy beam~^ for a beam 
size of 2'.'94x2'.'72; at 4.9 GHz, the noise was 15.0 /iJy beam"^ for a beam size of 5'.'09x4'.'43. 

A few radio sources with strengths of hundreds of microjansky to a few millijansky were 
found in the images, but we discuss only Gl in this Letterjf] At 8.4 GHz, an apparent source 
with a flux density of 28 ± 6 /iJy [corresponding to 2 x 10^^ W Hz~^ for distance modulus 
of (m — M) = 24.42 mag (Meylan et al. 2001)] was f ound approximately one arcsecond 



from the Gl optical position reported by iMeylan et al.l (120011 ): this radio source has J2000 
coordinates of a = 00''32'"46.54^ 6 = 39°34'39.2". Figure 1 shows our 8.4 GHz image of the 
20" by 20" region centered on Gl; thi s imag e includes a 1 sigma error circle of 1"5 radius 



for the X-ray position found by iKongI (120071 ). The radio position has an estimated error of 
0"6 in each dimension (not shown in the figure), derived by dividing the beam size by the 
signal-to-noise ratio0 

The a priori probability of finding a 4.5 a noise spike or background source so close 
to Gl is quite small, as indicated by the lack of any other contours of similar strength in 



^The image areas covered were much larger than the primary beams of the individual VLA antennas, in 
order to provide the best possible subtraction of confusing sources. 

■^The strongest nearby radio source is a 1 mJy (at 8.4 GHz) object located approximately 100" from Gl; 
this object and the other weak sources might be X-ray binaries or supernova remnants if located at the 
distance of M31, but their numbers also arc consistent with the possibility that some could be background 
extragalactic sources. 

^Transfer of the phase from the local phase calibrator makes an insignificant contribution to the radio 
source position error, relative to the error imposed by the limited signal strength. 



Figure 1. If we hypothesize that there are 9 independent beams (roughly 8" by 8") within 
which a source would be considered to be associated with Gl, then the probability of a 4.5 o 
noise point close to Gl is less than 10~^. Similar ly, the expected densit y of extragalactic 
radio sources at 28 /iJy or above is 0.25 arcmin"^ ( iWindhorst et al.lll993l ). or 4 x 10~^ in a 
box 8" on a side, making it unlikely that we have found an unrelated background source. In 
order to search for possible data errors that might cause a spurious source, we have subjected 
our data set to additional tests, imaging data from the two days separately, and also imaging 
the two different intermediate frequency channels separately. The Gl radio source remains 
in the images made from each data subset, with approximately the same flux density and 
position. The overall significance is reduced by 2^/^ to approximately 3 a in each image made 
with about half the data, as expected for a real source with uncontaminated data. Other 
2.5 cr-3 cr sources appear in the central 20" box in some subsets of half the data, consistent 
with noise statistics, but none is above the 3.5 a level in the full data set. Thus, all tests 
indicate that the detection of Gl is real, and we will proceed on that basis for the remainder 
of this paper. At 4.9 GHz, we find no detection at the Gl position, but the much higher 
noise level provides us only with very loose constraints on the source spectrum (see below). 



3. Origin of the Gl Radio Emission 



Merloni et al.l (120031 ) and lFalcke et al.l (120041 ) have quantified an empirical relation (or 



"fundamental plane" ) among X-ray and 5 GHz radio luminos i ty an d BH mass; we use the 
Merloni et aD J2003h relation Lr oc L^^M^^. iMerloni et aD hm'k s analyzed this relation 
in the context of accretion flows and jets associated with massive BHs. One might expect 
some general relationship among these three quantities, if an X-ray-emitting accretion flow 
onto a massive BH leads to creation of a synchrotron-emitting radio jet, with the detailed 
correlation providing some insight into the nature of that flow. By comparing the em pirically 
determined relation with expectations from theoretical models, iMerloni et al. J2003h deduced 
that the data for BHs emitting at only a few percent of the Eddington rate are consistent with 
radiatively inefficient accretion fiows and a synchrotron jet, but inconsistent with standard 
disk accretion models. 



Maccarond (120041 ) scaled the fundamental-plane relation to values appropriate for an 
IMBH in a Galactic globular cluster; we rescale their equation here to find a predicted radio 
fiux density of 
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Using the previously cited X-ray luminosity and BH mass for Gl and our adopted 
distance modulus, this predicts a 5 GHz flux density of 77 /iJy for Gl. However, taking into 
account the 30% uncertainty in the IMBH ma ss, the unknown sp ectral index of the radio 
emission, and the dispersion of 0.88 in logLpi (iMerloni et al.ll2003l ). the predicted 8.4 GHz 
flux density for Gl is in the range of tens to a few hundred microjansky. Thus, our radio 
detection of 28 /iJy at 8.4 GHz is consistent with the predictions for a 1.8 x 10^ Mq IMBH, 
but strongly inconsistent with a 10 Mq BH. Since neutron star X-ray bi naries in a variety 
of st ates have radio/X-ray ratios much lower than BH X-ray binaries (IMigliari fc Fender 
20061 ) . and thus another 2 orders of magnitude below the observed value, stellar-mass X-ray 
binaries of any type are ruled out as the possible origin of the radio emission in Gl. 

We can use the radio/X-ray ratio to assess other possible origins for the radio emission. 
Here, we use the ratio Rx = i'L,j{^A GHz)/Lx(2 — 10 keV) as a fiducial marker. For Gl, 
i?x ~ 5 X 10~^, which is considerably lower than Rx ~ 10~^ that is common to the Galactic 
supernova remnant Gas A, low-luminosity active galactic nuclei (supposing Gl might be a 
stripped dwarf elliptical galaxy), and most ultraluminous X-ray sources (cf. Table 2 of Neff, 
Ulvestad, & Gampion [2003], and references therein). 

It is of interest to compare the Gl source to various relatives of pulsars as well. For 
instance, Gl is within the wide ra nge of both luminosity and ra dio/X-ray ratio observed 
for pulsar wind nebula e (PWNs) (IFrail &: ScharringhausenI 119971 ). less luminous than the 
putative PWN in M81 jBietenholz. Bartel. fc Rupenl l 120041). but considerab ly more luminous 
than standard pulsars or anomalous X-ray pulsars (iHalpern et al.l l2005l ). The 8.4 GHz 
luminosity of Gl is similar to that of the magnetar SGR 1806 — 20 about 10 days after its 
outburst in late 2004, and the lack of a 4.9 GHz d etection would be con sistent with the fading 
of SGR 1806 — 20 two months after the outburst (IGaensler et al.ll2005l ). However, there is no 
published e vidence for a gamma-ray ou tburst from Gl, and the relatively steady apparent 
X-ray flux (jPooley fc RappaportI l2006l ) also argues against a transient source. Thus, the 
only stellar- mass object that might a ccount for the radio and X-ray emission would be a 
PWN; using the scaling law given by iFrail &: ScharringhausenI (119971 ). we find a likely size 
of ~ 10 milliarcseconds for a PWN radio source in Gl, implying that high- sensitivity VLBI 
observations could distinguish between a PWN and IMBH origin for the radio emission from 
Gl. 



Knowledge of the radio spectrum of Gl could provide more clues to the character of 
the radio emission, although either a PWN or an IMBH accretion flow might have a flat 
spectrum. In any case, our 5 GHz observation simply is not deep enough. If we choose 
a 2 cr upper limit of 30.0 /iJy at 4.9 GHz (2 a chosen since we know the position of the 
8.4 GHz source with high accuracy), we derive a spectral index limit of a > —0.12 ± 0.99 
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(for Su oc z/"*"", 1 cr error in spectral index), which has httle power to discriminate among 
models. 

The X-ray emission from Gl may b e due to Bondi accretion o n the IMBH, either from 
ambie nt clu ster gas or from stellar winds (jPooley fc Rappaportlbooeh . IHo. Terashima. fc Okaji 



ima 



(I2OO3I ) and iPooley fc RappaportI (120061 ) give approximate relations for the Bondi accretion 
on an IMBH in a globular cluster; for an ambient density of 0.1 cm~^, an ambient speed of 
15 km s~^ for the gas particles relative to the IMBH, and a radiative efficiency of 10%, the 
Bondi accretion luminosity for the Gl I MBH would be ~ 3 x 10^^ erg s s~^. The X-ray lumi- 
nosity of 2 X 10^^ ergs s~^ measured by iPooley fc RappaportI (120061 ) thus implies accretion 
at just under 1% of the Bondi rate. Given that Lx/L^dd ~ 10~^, a more likely scenario is 
that Gl accretes at closer to 10% of the Bondi rate but with a radiative efficiency under 1%. 
In this context, we note that the radio/X-ray ratio for Gl is logi?x > —4.3, which is above 
the v alue of —4.5 used to divide radio-quiet from radio-loud objects (ITerashima fc Wilson 
2003h FI Gl therefore should be considered radio-loud, as inferred for BHs in galactic nuclei 
that radiate well below their Eddington luminosities (Ho 2002). 

If the globular clusters in our own Galaxy also have central BHs that are 0.1% of their 
total masses, and they accrete and radiate in the same way as Gl, many would have expected 
5 GHz radio flux densities in the 20-100 /xJy range; flux densi ties often would be in t he 1- 
10 jUJy range ev e n for less efficient accretion and radiation (jMaccarone et al" 2005). As 
Maccarone et al.l (120051 ) summarize, t here are few r adio images of globular clusters that go 
deep enough to test this possibility. iFenderl (120041 ) points out that the Square Kilometer 
Array (SKA) will be able to test for the existence of IMBHs in many globular clusters. 
However, based on our results for Gl, we suggest that it is not necessary to wait for the 
SKA; t he current VLA can reach the hypothesized flux densities with some effort. The 
EVLA (jUlvestad et al.ll2006l ). scheduled to be on line in about 2010, will have 40 times the 
bandwidth and 6.3 times the sensitivity of the current VLA in the frequency range near 
8 GHz. This will enable the EVLA to reach the 1 yuJy noise level in approximately 12 hours 
of integration, thus probing the range of radio emission predicted by lMaccarone et al.l (120051 ) 
for many globular clusters. 



lower limit is given for i?x because this quantity traditionally is given in terms of the 2-10 keV 
luminosi ty, whereas log-Rx = ~4-3 w ould correspond to the value computed for the 0.2-10 keV luminosity 
given bv lPoolev &: RappaportI (|2006h . 
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4. Summary 

We have detected faint radio emission from the object Gl, a globular cluster or stripped 
dwarf elliptical galaxy in M31. The emission has an 8.4 GHz power of 2 x 10^^ W Hz~^. As- 



sumi ng that the radio source is associated with the X-ray source in Gl ( iPooley &: Rappaport 



20061 ). the radio/X-ray ratio is consistent with the value expected for an accreting ~ 2 x 
10^ Mq BH. Thus, the radio detection lends support to the presence of such an IMBH within 
Gl. The other possible explanation, a pulsar wind nebula, could be tested by making very 
high-sensitivity VLBI observations of Gl. 
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Fig. 1. — VLA C configuration 8.4 GHz image of tlie vicinity of Gl in M31, with contours 
at intervals of 0.5 times the r.m.s noise. The lowest contour is at 3 times the noise of 
6.2 /iJy beam~^, and negative contours are shown dashed; the synthesized beam is shown in 
the box in the lower-left corner. The (0,0) point is at a J2000 position of a = 00'*32™46.60^, 
6 = 39°34'40.0". The l?5 -radius circ le about this central point represents the 1 a error in 
the X-ray position of Gl (IKongl 120071 ). 



